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Abstract. The van der Waals and Casimir interactions of a hydrogen atom 
(molecule) with a single-walled and a multiwalled carbon nanotubes are compared. 
CZ3 . It is shown that the macroscopic concept of graphite dielectric permittivity is already 

applicable for nanotubes with only two or three walls. The absorption of hydrogen 
atoms by a nanotube at separations below one nanometer is considered. The lateral 
force due to exchange repulsion moves the atom to a position above the cell center, 
where it is absorbed by the nanotube because the repulsive force cannot balance the 
^ • van der Waals attraction. 
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PACS numbers: 73.22.-f, 34.50.Dy, 12.20.Ds 
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1. Introduction 

! Considerable recent attention has been focused on applications of the van der Waals and 
Casimir forces in nanotechnology. When the characteristic sizes of microdevices shrink 

00 . 

O , below a micrometer, the collective quantum phenomena caused by the existence of zero- 
point oscillations of the electromagnetic field come into play. At separations below 
100 nm the role of physical phenomena originating from vacuum oscillations (primarily 
of the van der Waals and Casimir forces) can match the role of characteristic electric 
force and may even become dominant. This was first stressed long ago in pQ and 
was commonly accepted in the beginning of the Third Millenium. By then several 
experiments on measurement of the Casimir force between metallic surfaces had already 
been performed [2]. This laid the groundwork for the application of dispersion forces 
(which is a generic name for the van der Waals and Casimir force) in nanomechanical 
devices [31 IH [5], noncontact atomic friction P, [7J EJ [91 HO], carbon nanostructures 
[T2| [T3l [H] and related subjects. Thereafter, the precision of measurements of the 
Casimir force was significantly increased [15] and different methods to control the force 
magnitude were elaborated [33] opening possible applications to nanotweezers, nanoscale 
actuators and other nanomachines. 
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Dispersion interaction is also of much importance for the understanding of 
absorption phenomena of atoms and molecules by nanostructures, This subject is 
currently topical in connection with the problem of hydrogen storage in carbon 
nanotubes. Until recently it was studied using the phenomenological density-functional 
theory (see, e.g., [TU [121 E]). In [Tf[ [18] the Lifshitz theory of the van der Waals and 
Casimir force [19] was applied for the cases of an atom (molecule) interacting with a 
plane surface of a graphite plate or multiwall carbon nanotube. This was achieved by 
using a classical model of the frequency dependent dielectric permittivities of a uniaxial 
crystal and proximity force approximation [20] . 

Some properties of the monoatomic sheet of C atoms (graphene) admit a simplified 
model description in terms of the two-dimensional free electron gas. In doing so 
the graphene sheet is characterized by some typical wave number K determined by 
the parameters of the hexagonal structure of graphite. In [21] the interaction of 
the electromagnetic oscillations with such a sheet was considered and the reflection 
coefficients were found. In [22] the Lifshitz-type formula was obtained for the van der 
Waals and Casimir interaction between the two parallel plasma sheets. Using this model 
in [23] the interaction between graphene and a material plate, graphene and an atom 
or a molecule, and between a single-walled carbon nanotube and a material plate was 
also described by means of the Lifshitz-type formulas. Finally, in [23] the Lifshitz- 
type formula was obtained for the van der Waals and Casimir interaction between an 
atom (molecule) and a single-walled carbon nanotube. All above mentioned results are 
applicable if the separation distances are sufficiently large (typically larger than lnm). 

In this paper, we compare the van der Waals interactions of a hydrogen atom 
(molecule) with a single-walled and a multiwalled carbon nanotube. The former is 
approximately modeled as a cylindrical plasma sheet, whereas the latter as a cylindrical 
shell of finite thickness characterized by graphite dielectric permittivities for the ordinary 
and extraordinary rays. We arrive at the conclusion that the macroscopic concept 
of dielectric permittivity is already applicable for nanotubes containing only two or 
three walls (Section 2). We also consider the interaction of atoms or molecules with 
carbon nanotubes at separations below 1 nm where the Lifshitz-type formulas obtained 
in [2~2| [23] [24] are not applicable. Using the method of phenomenological potentials 
and disregarding the role of chemical forces, we find that at separations below 1 nm the 
exchange repulsion gives rise to the lateral force that moves hydrogen atoms towards 
the cell center. In the position above a cell center, the repulsive force cannot balance 
the van der Waals attraction. As a result, the atom penetrates inside the nanotube 
(Section 3). This effect is analogous to the breaks of constant force surfaces that arise 
when scanning the monoatomic tip of an atomic force microscope above a closely packed 
lattice in contact mode [25j [26]. The discussion of the obtained results and conclusions 
are contained in Section 4. 
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2. Comparison of atom-nanotube interaction in the cases of multiwalled 
and single-walled nanotubes 

The free energy and force between an atom (molecule) and a multiwalled or single-walled 
carbon nanotube separated by a distance a can be represented in the form 

~~/ m % Cs(a,T) m . Cp(a,T) . , 

F{a, T) = F(a, T) = (1) 

Here, T is the temperature of graphite which is supposed to be in thermal equilibrium 
with the surroundings. The coefficients C%{a,T) and Cf(o,,T) are defined in such a 
way that at short separations equation (pQ) reproduces the nonrelativistic van der Waals 
interaction (in this limit C3 and Cp do not depend on separation and temperature and 
it holds Cp = 3C3). The Lifshitz-type formula for the coefficient C^{a, T) was obtained 
in [TTJ, [23] by using the proximity force approximation [201 12Z] 



k B T R 4B + 3a 
C ^ T) = —yR^a\2jRTaj a{0) (2) 
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2r TM (iCi, y) + -^(r TE (iO, I/) - ^tm(iC/, 1/) 

Here, i? is the nanotube radius, a(u) is the dynamic polarizability of an atom or a 
molecule, ks is the Boltzmann constant and the dimensionless Matsubara frequencies 
(1 are connected with the dimensional ones by the equalities 

= -, 6 = 2tt^Z, ^ = t^, Z = 0,1,2,.... (3) 

In our model description, the reflection coefficients for the two independent polarizations 
of the electromagnetic field, ttm and ?"te; have different forms in the cases of multiwalled 
and single-walled carbon nanotubes. For multiwall nanotubes, they are expressed in 
terms of two dissimilar dielectric permittivities of graphite, e x (u) = e y (u>) and e z (u), 
where the crystal optical axis z is perpendicular to a cylindrical surface 

e x ie z iy 2 - ft 



rTu(Ki,y) 



e x ie z iy 2 + P z + 2 v /e xl e zl yf z coth[f z d/(2a)} 



fx - y 2 

r ™ V) = P x +y 2 + 2yf xC otWxd/(2a)\ (4) 

Here, the thickness of a nanotube is d = 3.4(iV — 1) A, where iV is the number of walls, 
and the following notations are introduced: 

£xi = e x (iQu c ), e z i = e z (i(iuj c ), (5) 
/, 2 = Z/ 2 + C 2 (^-l), f 2 x =y 2 + (?(e xl -l). 

The dielectric permittivities of graphite along the imaginary frequency axis are 
calculated in [17] on the basis of tabulated optical data from [28] . 
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For single-walled nanotubes, we use the simplified model description by means of 
a cylindrical plasma sheet. Then the reflection coefficients are expressed in terms of 
the wave number of the sheet, K, which is determined by the density n of ir electrons. 
Bearing in mind that there are two ir electrons per one hexagonal cell with a side length 
ro = 1.421 A, we arrive at 

K = 2tt — - = 6.75 x 10 5 m -1 , n = — ■=—. (6) 
mc 2 Zy/Zrl 

where e and m are the electron charge and mass, respectively. The resulting coefficients 

are 

2yaK 2aK 
2yaK + C?' ™ [Kh V) ~ 2^K + y 
The Lifshitz-type formula for the coefficient Cp(a,T) in §T§ can be presented using 
the same notations [171 [21]. We emphasize that within the considered models equation 
(0Q) with the coefficients C3 and Cp is applicable at both short and large separations, i.e., 
both in nonrelativistic and relativistic regimes, and also in the transition region between 
the two regimes. We, however, consider a > 1 nm because at shorter separations the 
atomic structure of nanotube wall and also other forces of different physical nature in 
addition to dispersion interaction should be taken into account. In the framework of used 
models, the error introduced from the application of the proximity force approximation 
is less that 1% within the separation region from to R/2 [27]. It is notable also that 
the above models do not take into account nanotube chirality. This can be included 
by using the optical data for nanotube complex index of refraction. We would like to 
underline that the used models, especially the model of a cylindrical plasma sheet in 
application to single- walled carbon nanotubes, do not claim a complete description of all 
nanotubes properties. Specifically, it remains unclear if it possible to describe nanotubes 
with surfaces like metals and like semiconductors by varying only one parameter K in 
the reflection coefficients ([7]). 

To perform computations using above equations in the case of hydrogen atoms and 
molecules, one neads the explicit expressions for the atomic and molecular dynamic 
polarizabilities. As was shown in [TTJEE], the dynamic polarizability can be represented 
with sufficient precision using the single-oscillator model 

= c2 , oi m {i£{) = 9m (8) 

+ it + it 

for hydrogen atom and molecule, respectively. Here, a a (0) = 4.50 a. u., u a = 11.65 eV, 
ct m (0) = 5.439 a.u., u a = 14.09 eV. 

We have performed computations of the van der Waals free energy of hydrogen 
atoms and molecules and multiwalled and single-walled carbon nanotubes modeled as 
discussed above. For multiwalled nanotubes, equations ([I]), (j2J), d3J, (|EJ) were used. For 
single- walled nanotubes equation was replaced with ([7j). The computational results 
for the coefficient C3 as a function of the number of walls are presented as solid dots 
in figure 1(a) (hydrogen atom) and in figure 1(b) (hydrogen molecule) for nanotubes 
with R = 5 nm. In both cases the solid dots connected with a solid line are related 
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Figure 1. The van der Waals coefficient as a function of the number of walls for 
hydrogen atom (a) and molecule (b) interacting with the multiwalled carbon nanotube 
(solid dots connected with solid lines) and with the single-walled carbon nanotube 
(solid dots 1, 2 and 3) both of R — 5 nm radius spaced at 1, 2, and 3 nm from the atom 
(a) or molecule (b). The thicknesses of multiwalled nanotubes with 2, 3, 4 and 5 walls 
are 0.34, 0.68, 1.02 and 1.36 nm, respectively. 



to a multiwalled nanotube and the solid dots marked by 1, 2, and 3 are related to a 
single-walled nanotube for separation distances of a = 1, 2, and 3nm, respectively. For 
a "multiwalled" nanotube with only one wall C3 = 0, as expected, because in this case 
nanotube thickness d = and hence the reflection coefficients (J4]) vanish. As is seen in 
figure l(a),(b), the magnitudes of the van der Waals coefficients C3 for a multiwalled 
nanotube with two walls at a separations of 2 or 3nm from an atom or a molecule are 
in correct relation to the magnitudes computed for a single-walled nanotube (solid dots 
marked 2 and 3). Thus, the macroscopic concept of graphite dielectric permittivity is 
already applicable for nanotubes with only two or three walls provided the separation 
distance to an atom or a molecule is larger than 2 nm. For nanotubes with three walls 
the same holds for sepataions of 1 nm between an atom or a molecule and a nanotube. 
These results are expected because it is known that surface corrections to local fields in 
ordinary crystals generally become negligible already in the third or fourth lattice plane. 
Similar results but with other numerical values are also valid for the force coefficient 
Cf- To give an example, for nanotubes of R = 5 nm radius at a = 1 nm from a hydrogen 
atom C3 = 0.0585 a.u. but Cp = 0.197 a. u. At the same separation but for nanotubes 
of R = 2nm radius, C3 = 0.0503 a.u. and Cf = 0.175 a.u. 

3. Absorption of hydrogen atoms by carbon nanotubes 

The above formalism based on the Lifshitz-type formulas permits to calculate the 
attractive van der Waals and Casimir force acting between both multiwalled and single- 
walled carbon nanotubes and hydrogen atoms or molecules down to lnm separation 
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distance using the simplified models described above. However, at shorter separations 
the atomic structure of nanotube and some other physical interactions in addition to 
dispersion forces should be taken into account. Because of this, complete calculation 
of absorption power of carbon nanotubes requires a detailed investigation of the 
microscopic interaction mechanisms including chemical forces and short-range exchange 
forces. The comparative role of different forces in absorption process has not been yet 
investigated. Below we consider the role of dispersion and exchange forces at separations 
less than lnm. We argue that the combined action of the interatomic attractive van 
der Waals force and of the repulsive exchange force, taken into account using a simple 
model, leads to the absorption of hydrogen atom by a nanotube. The investigation of 
the role of chemical forces is an interesting problem to be solved in future. 

The exchange repulsion between C and H atoms at a separation r apart can be 
approximately described by means of the phenomenological potential U Tep (r) = a/r 12 
|30j, where a is some coefficient. Papers [251 I2S] investigated the scanning of the 
monoatomic tip of an atomic force microscope along the constant force surface above 
the closely packed crystal lattice in contact mode. It was supposed that the repulsive 
potential U rcp (r) is the single factor to be taken into account. As was shown in 
[25| [26] , if the initial height of the tip h satisfies the condition h/r > 0.61 (here, 
r is the equilibrium interatomic distance in the lattice), the constant force surface 
under consideration is continuous. If, however, the initial height is sufficiently small 
(h/r < 0.61), the constant force surfaces have breaks above the regions between crystal 
atoms. For graphite lattice the atom in the center of a hexagonal cell is missing which 
only increases the areas of the breaks. 

Let us consider now the hexagonal lattice shown in figure 2(a) with a hydrogen atom 
at a height h above it in the positions between points A and B (the lattice parameter is 
r ). There is the lateral force acting on an atom at a height h directed to a cell center. 
This can be seen in the following way. Let us consider the total repulsive energy due 
to interaction with a potential U vep with all nearest C atoms with coordinates (xj, t/i, 0), 
e.g., with 4 atoms in the position A and 6 atoms in the position B (the first coordinative 
sphere). This energy is given by 



where x = in the position A and x = r in the position B. In figure 2(b) we plot the 
ratio E iep (x, h)/E Tep (0, h) as a function of (r — x)/r Q for different heights: h = 1.5r 
(line 1), h = ro (line 2) and h = 0.5ro (line 3). As is seen in figure 2(b), the equilibrium 
position of an atom at x = (point A) is unstable, whereas the equilibrium position 
of an atom at x = ro (point B) is stable. Thus, atom will be displaced by a horizontal 
force to the position B above the cell center (in fact it should be displaced to the center 
of one of the three hexagonal cells around a position A). 

We further suppose that the total interaction of H atom with C atom at a separation 




(9) 
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Figure 2. (a) Hydrogen atoms in the positions A and B at a height h above the 
hexagonal lattice of C atoms, (b) The relative energy of exchange repulsive interaction 
between the hydrogen atom at a height h above the hexagonal lattice of C atoms and 
neighboring C atoms as a function of (ro — x)/tq. Here x is the lateral displacement 
from the C atom below a point A to the center of the cell. 



r below 1 nm is described by the Lennard- Jones interaction potential 

UUr) = £ " (10) 

where (3 is the constant of interatomic van der Waals interaction. The equilibrium 
separation distance h$ between one C and one H atom satisfies the condition 
dU LJ (r) Ua 6(3 Q 



dr r =h h o 3 h l 
This leads to the relation (3 = 2a/ h®. 

The total force acting on an H atom at a height h in the position B is determined 
by the six neighboring C atoms. It is equal to the sum of vertical projections of forces 
with potential (TIDi) 

72ah 36(3h 
B {rl + Wy (r 2 + ^)4- [ > 

The hydrogen atom would be absorbed by a nanotube if Fb < for any value of h. 
Taking into account the above connection between a and (3, we obtain from (fl2l) 

i < 0. (13) 



(r 2 + /i 2 )3 hi 
This inequality is for sure satisfied if h < r . 

It is common knowledge [301 EI] that the equilibrium position of the two atoms 

h Q ^Rfi + Rfl<R^ + B%\ (14) 

where R^ n and R^ are the so-called ionic and atomic radii, respectively. For C and H 
atoms under consideration it holds R^ = 0.77 A and R^ = 0.53 A[32]. Thus, in our 

case 

h <RP + RW = 1.3 A. (15) 

If to take into account that for graphite cell r = 1.42 A, we arrive at the conclusion that 
the condition ho < ro is satishied and hydrogen atom will be absorbed by a nanotube. 
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4. Conclusions and discussion 

In the above we have presented the Lifshitz-type formulas for the van der Waals and 
Casimir free energy and force for the configuration of hydrogen atoms or molecules 
in close proximity to multiwalled or single-walled carbon nanotubes. Multiwalled 
nanotubes are approximately described by graphite dielectric permittivity, whereas 
single-walled nanotubes are considered in the approximation of a two-dimensional gas 
of free electrons. Both descriptions are of model character and do not claim complete 
description of all nanotube properties. They are applicable at separations larger than 
1 nm. By comparing of calculation results for multiwalled nanotubes with those for 
single-walled nanotubes, the conclusion was made that the macroscopic description 
using the concept of dielectric permittivity is already applicable for nanotubes with 
only two or three walls depending on the separation distance to the hydrogen atom. 

We have also considered separations below 1 nm where the exchange repulsion plays 
an important role. By disregarding the role of chemical forces, it was shown that under 
the influence of a lateral force originating from exchange repulsion, the hydrogen atom 
is moved towards a cell center. Simple analysis shows that at atomic distances above 
the cell center the exchange repulsion cannot balance the van der Waals attraction. As 
a result the hydrogen atom is absorbed by the nanotube. This is analogous to the effect 
of breaks on the constant force surfaces which arise when scanning of the monoatomic 
tip of an atomic force microscope above a surface in the strong repulsive mode. In future 
it would be interesting to consider the influence of the previously absorbed atoms on 
the absorption process and to determine the resulting absorption rate. 
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